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Abstract 

q . The equation of state for quark-gluon plasma is obtained, being valid for arbitrary 

values of the coupling parameter and of temperature. This equation is constructed 
from perturbative expansions for free energy and for renormalization function. Sum- 
mation of asymptotic perturbative series is accomplished by means of the self-similar 
approximation theory. 
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1 Introduction 

o. 

Under an equation of state one implies the dependence of pressure or free energy on 
Ph| some characteristic parameters, for instance, coupling parameter or temperature. The 

Q_i' definition of pressure is 

p=-ln Tr e~ / , 



X 



where T is temperature; V, volume; and H is the Hamiltonian defining the considered 
system. In recent years there has been much attention paid to understanding the equation 
of state in quantum chromodynamics. The state of matter, called quark-gluon plasma, 
almost certainly existed in the early Universe up to 10 _5 s after the Big Bang and is likely 
to be found in the interior of neutron stars. One hopes that this state can be formed 
in the collisions of relativistic heavy nuclei where it can be studied in a controlled and 
systematic way [1-6]. The formation of quark droplets in colliding nuclei can occur at 
already existing accelerator energies, leading to the cumulative effect [7]. 

There are two ways of obtaining the equation of state in quantum chromodynamics. 
One is a numerical way based on lattice simulation, and another way is by employing 
statistical models. A discussion of these ways can be found in the recent review [6]. In 
this report we advance the third possibility of deriving the equation of state in quantum 
chromodynamics. This novel way is based on the summation of perturbative series in 
powers of the coupling parameter. The method of summation we use here is a variant of 
the self-similar approximation theory [8-15] employing self-similar exponential approxi- 
mants [16]. This method has been successfully used for constructing accurate equations 
of state for several physical systems [15,16]. 
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2 Summation of series 



Perturbation theory can be used in QCD at high temperatures when the coupling pa- 
rameter becomes small. One usually considers the case of zero chemical potential and 
assumes that the temperature is high enough that fermion masses can be ignored. An 
expansion for the QCD free energy or pressure, in powers of the coupling parameter g has 
been known to fourth order [17-20] and recently was calculated to fifth order [21,22]. 

The perturbative expansion for the pressure of high-temperature gauge theory with 
massless fermions in four dimensions can be written [21] as 



p(g) 



w 2 dj 



T 4 £(c n + c' n In g)g», 



(1) 



where the coefficient are 
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Here the dimensionless regularization is used, the scale \x corresponds to the modified min- 
imal subtraction scheme, ((■) is the Riemann zeta function, 7^ is the Euler-Mascheroni 
constant, and for SU (N c ) theory, with N c colours and with rif fermions in the fundamental 
representation, 

d A — N 2 — 1 , d F = N c n f , C A = N c , 



1 



N 2 



>1F 



AN a 



For QCD with A^ c = 3, one has 
d A = 8 dp — 3n/ , 



C A = 3 , 
2 



n f . 



S- 



2F 



Then expansion (1) can be written as 

p(g) ~ ^T 4 c n «f + In <7 £ cU") , (2) 
with the nonzero coefficients 

3/2 



c = 1 + ^n/ , c 2 = -0.09499 (l + ^n/) , c 3 = 0.12094 (l + i 



/i 
T 



c 4 = 0.04331 (l + j^n/) In (l + j^n/) + 
+0.01733 - 0.0076372/ - 0.00088nJ - 0.01323 + ^n/) (l - — n/) In 
c 5 = - f 1 + -n/J (0.12806 + 0.00717n/ - 0.00027nj) + 

+0.02527 (l + i % ) 3/2 (l - A n/ ) In | , 

c 4 = 0.08662 + in/ 
It is convenient to introduce the dimensionless function 

*)^( 1 + l-) r4 ' < 3 » 

being pressure (2) reduced to the Stefan-Boltzmann limit p(0). As g — > 0, the reduced 
pressure (3) reads 

/G?)^E^" + lns£<^, (4) 

n=0 ra=2 

where a n = c„/c and = c' n /c . Truncating the asymptotic series (4), we have the 
sequence of perturbative approximations 

f k (g) = 1 + a 2 g 2 + ... + a k g k , (5) 

for k — 2, 3, 4, 5, and 

U (9) = 1 + a 2 <? 2 + • • • + « 6 <? 6 + In 5 (4<7 4 + a^ 6 ) (6) 

for k = 6. Note that the term g 4 \ng has to be ascribed to the sixth-order approximation 
since only then one gets the series factoring Ing, which could be summed. 

The renormalization scale u can be considered as a control function which can be 
defined from a fixed-point condition, e.g. from the minimal-difference condition [23] 

fk(9,fJ>) ~ fk-i(g,V>) = , (7) 

applied to the first k approximation when the scale g appears. In our case, this is k — 4, 
for which condition (7) is equivalent to 04 = or C4 = 0. Then we find 

a 0.04331 (l + ±n/) In (l + ±n/) + 0.01733 - 0.00763n/ - 0.00088n^ 

In — = In 7 = 7 — r— 7 r . (8) 

0-01323 (l + ±n f ) (l - |n/) 
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This gives /i ~ T, which is physically reasonable since T is the natural typical scale for 
the quark-gluon plasma. For example, 7 = 3.70580 if nf = 0; 7 = 1.69564 is rif = 3, 
and 7 = 0.99696 if rif — 6. Note that if one would take = 5 in condition (7), that is 
a 4 = c 4 = 0, one would get a physically unreasonable value ji ~ 100T. 

Following the method of self-similar exponential approximants [15,16], starting from 
the sequence {fk(g)} with fk(g) given in (5), we construct the exponentials 

F 2 (g, r) = exp (b 2 g 2 r) , 

F 3 (g, r) = exp (b 2 g 2 exp (b 3 gr)) , (9) 

F 5 (g, t) = exp (b 2 g 2 exp (b 3 g exp (& 5 # 2 r))) , 

in which F 4 = F 3 because of (7), the coefficients are 

, c 2 a 3 c 3 a 5 c 5 

6 2 = 02 = — , o 3 = — = — , 05 = — = — , 
c a 2 c 2 a 3 c 3 

and t = Tk(g) is a control function for the corresponding expression F k (g, r), so that 

F k (g,r) - F fe _!(5,r) = , T = r k (g) , 

similarly to (7), with T 2 (g) = 1. 

After defining the control functions r k (g), we obtain the self-similar approximants 

ftig) = F k (g,r k (9)) • (10) 

In particular, 

f 2 (g) = exp(b 2 g 2 ) , 
f* 3 (g) = ex P (b 2 g 2 r 3 (g)) , (11) 

/a (#) = ex P ( b 2g 2 exp (b 3 gT 5 (g))^j , 
with the control functions defined by the equations 

r 3 = exp(fc 3 #T 3 ) , r 5 = exp(& 5 si 2 T 5 ) . 

In expansion (4) and, respectively, in approximants (11), the coupling parameter g = 
g([i) is the running coupling described by the renormalization group equation 

^ = (3(g) ■ (12) 

The renormalization function (3(g), as g — > 0, has the expansion [24] 

(3(g) ~ ^ 3 + /% 5 + /V, (13) 



in which 
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As the self-similar approximants for (13), we have 



(3*(g) = ag 3 exp(bg 2 ) , 

fi(g) = ag 3 ex V (bg\ 7 (g)) , (14) 

where the control function r 7 {g) is the solution of the equation r 7 = exp(cg 2 r 7 ) and the 
coefficients are a = /3 3 , b = P5/P3, and c = /3 7 //3 5 . 

To solve equation (12) with the right-hand side given by one of the functions from 
(14), we need an initial condition. For the latter, we may accept the value a s {mz) = 0.119 
of the strong coupling constant a s = g 2 /Atx at the Z° boson mass m z [25]. Thus, the 
initial condition for equation (12) is 

g(m z ) = 1.222856 , m z = 91.187 GeV . 

Solving (12), we find g = g(fi). Substituting the latter in (10) and using the relation 
H = 7T following from (8), we obtain the reduced pressure 

fk(T) = f*MlT)) (15) 

as a function of temperature. 

3 Results and Conclusion 

The results of our calculations for rif = 6 are presented in Figs. 1 to 3. The self-similar 
approximants /£ (g) given by (11) for the reduced pressure (3), as functions of the coupling 
parameter g, are shown in Fig. 1. The solution g(fi) to the evolution equation (12) is 
drawn in Fig. 2. Both functions in (14) give practically the same solution g(/i). Figure 3 
presents the equation of state (15). 
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Figure 1: The behaviour of the reduced pressure f(g) = p(g)/p(0) presented by the 
self-similar approximations /2(g) (solid line), f£(g) (long-dashed line), and f$(g) (short- 
dashed line). 
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Figure 2: The running coupling g(/i) as a function of the scale /i, obtained from the 
renormalization-group equation (12). 



f(T) 




Figure 3: The reduced pressure f(T) given by the self-similar approximants f2(T) (solid 
line), fz(T) (long-dashed line), and fe{T) (short-dashed line). 
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The overall behaviour of the reduced pressure f(T) is in agreement with that found 
in lattice simulations and in statistical modelling (see review [6]). At the temperature 
T c pa (150 — 200)MeV the pressure sharply drops down to practically zero, which can 
be interpreted as confinement. The curves corresponding to subsequent approximations 
k — 2, 3, 5 are close to each other. The qualitative behaviour of f(T) for different numbers 
of flavours in the diapason < rif < 8 is the same, with increasing T c as n/ decreases. 
Certainly, in the vicinity of T c ~ 200MeV one should take into account the quark masses, 
especially those of heavy quarks. But what is very interesting is that, starting from pertur- 
bative expansions in the high-temperature region of asymptotic freedom, we obtained an 
equation of state with the qualitatively correct behaviour, at all temperatures, including 
the existence of confinement at T c ~ 200MeV. 
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